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Abstract. We present results from a kinematical study of the gas in the nucleus of a sample of three LINER galaxies, obtained 
from archival HST/STIS long-slit spectra. We found that, while for the elliptical galaxy NGC 5077, the observed velocity 
curves are consistent with gas in regular rotation around the galaxy's center, this is not the case for the two remaining objects. 
By modeling the surface brightness distribution and rotation curve from the emission lines in NGC 5077, we found that the 
observed kinematics of the circumnuclear gas can be accurately reproduced by adding to the stellar mass component a black 
hole mass of Mbh = 6.8I28 ^ 10* M© (uncertainties at a Icr level); the radius of its sphere of influence (iJjph ~ 0'.'34) is well- 
resolved at the HST resolution. The BH mass estimate in NGC 5077 is in fairly good agreement with both the Mbh - (with 
an upward scatter of ~ 0.4 dex) and Mbh - cr correlations (with an upward scatter of 0.5 dex in the Tremaine et al. form and 
essentially no scatter using the Ferrarese et al. form) and provides further support for the presence of a connection between the 
residuals from the Mbh -o" correlation and the bulge effective radius. This indicates the presence of a black hole's "fundamental 
plane" in the sense that a combination of at least cr and drives the correlations between Mbh and host bulge properties. 
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1. Introduction 

It is now clear that the presence of a supermassive black 
hole (SMBH) is a common, if not universal, feature in 
the center of galaxies. In fact, since active galactic nuclei 
(AGNs) are thought to be powered by mass accretion onto 



and bulge velocity dispersion crh„i ([Ferrarese & Merritt I2OOOI: 
Gebhardt et al]|2000l: iTremaine et allbOoHl 



tivity in nearby 


galaxies (Heckman 1980'; 'Maoz et al."l995'; 


HoetaL 1997allb 


; Braatz et al...l997; Barth et al...l998. .1999; 


Naearetal. 2002) has led to the conclusion that a signif- 



icant fraction of galaxies in the Local Universe must host 
SMBH. This conclusion is now supported by direct measure- 
ments of SMBH masses in the center s of nearby galaxies ob- 



tained with different techniques (see iFerrarese & FordI 12005 
for a review). These measurements indicate that the BH mass 
Mbh is related to the properties of t he host galaxy, such as 
bulge luininosity Lhui and ma ss Mbui ( Kormendv & Ric hstone 
I995I; iMagorrian et aP 1 19981: iMarconi & Hunt 2003.) ligh t 



1995t IMagorrian et al.1 119981: IMarconi & Hunt 2003.) ligh 
concentration Joraham et al.l bOOll iGraham & Driveil l2007h 
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The existence of any correlations between Mbh and host 
bulge properties supports the idea that the growth of SMBHs 
and the formation of bulges are clo sely linked (Silk & ReesI 
19981; iHaehnelt & Kauffmannllioool) . This has profound im- 
plications for the process of galaxy formation and evolution. 
Moreover, SMBH mass estimates inferred via the above corre- 
lations, when more direct methods are not feasible, enter into 
a variety of important studies spanning from AGNs physics to 
the coeval formation and evolution of the host galaxy and its 
nuclear black hole. 

These results need, however, to be further investigated by 
increasing the number of accurate BH mass determinations in 
nearby galactic nuclei to set these correlations on a stronger 
statistical basis. In particular, such a study has the potential of 
establishing the precise role of the various host galaxy's param- 
eters in setting the resulting BH mass. To date, reliable SMBH 
detect ions have been obtained for a limited number of galaxies 
(~ 30. Ferrarese & Ford.2005 ), with the bulk of Mbh estimates 
in the range of 10^ - IO'Mq. To add reliable new points to 
the Mbh- host galaxy's properties planes is then a fundamen- 
tal task for future developments of astronomical and physical 
studies. 
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One widely applicable and relatively si mple method of de 
tecting BHs is based on gas kinematics (e.g. Harms et"al]|l994 



iFerrarese et al.lll996HMacchetto et al.lll997tlBarth et al.ll200lb . 
through studies of emission lines from circumnuclear gas disks, 
provided that the gas velocity field is not significantly influ- 
enced by non gravitational motions. However, the purely grav- 
itational kinematics of the gas can be established a posteriori 
from successful modeling of the gas velocity field under the 
sole influence of the stellar and black hole potential. 

The Space Telescope Imaging Spectrograph (STIS) on- 
board HST is still the most suitable instrument for such studies 
as it provides a high angular resolution (~ (X'l) in the optical 
spectral region. In this band brighter emission lines are found 
with respect to the infrared, the only band accessible at high 
resolution with ground-based adaptive optics telescopes. The 
wealth of unpublished data contained in STIS archives repre- 
sents an extraordinary and still unexplored resource. With the 
aim of finding galaxy candidates to provide a successful SMBH 
mass measurement, we performed a systematic search for un- 
published data in the HST STIS archive. 

In this paper we present the results obtained for a sam- 
ple of three LINER galaxies. The galaxies were observed with 
STIS on HST during Cycle 7 under Proposal ID 7354 and were 
part of a larger sample of eight selected obj ects. The results on 
NGC 3998 have already been published dOe Francesco et aT 



20061 hereafter Paper I). In a forthcoming paper, we will com- 



plete the proposal targets and discuss the relation between BH 
mass, host galaxy properties, and nuclear activity. Table [T]lists 
the galaxies and their main physical properties. Basic data are 
from the Lyon/Meudon Extragalactic Database (HyperLeda) or 
NASA/IPAC Exti-agalactic Database (NED). Distances are cal- 
culated with Ho = 75 km s"' Mpc"' and are corrected for Local 
Group infall onto Virgo. 

The paper is organized as follows. In Sec. |2] we present 
HST/STIS data and the reduction that lead to the results de- 
scribed in Sec.|3] In Sec.|4]we model the observed emission- 
line rotation curve for NGC 5077 and show that the dynam- 
ics of the circumnuclear gas can be accurately reproduced by 
circular motions in a thin disk when a point-like dark mass is 
added to the stellar potential. Our results are discussed in Sec. 
|5] and summarized in Sec.|6] 

2. HST data and reduction 

The three galaxies considered were observed with STIS on 
HST with the G750M grating and the 52"x O'.'l sht. Data were 
acquired at different slit positions, following a perpendicular- 
to-slit pattern with a step of U.' 1 and the central slit centered on 
the nucleus. The spectra obtained for each source, NUC for the 
nuclear slit, Nl-Nn (to North), and Sl-Sn (to South) for the off- 
nuclears, were retrieved from the public archive. Table [T] lists 
the number of slit positions, the orientation of the slit (from 
north to east), the exposure time for each positioning, and the 
observation date for the galaxies of the sample. 

The data were obtained with a 2x1 on-chip binning of 
the detector pixels and automatically processed through the 
standard CALSTIS pipeline to perform the steps of bias and 
dark subtraction, applying the flat field and combining the two 



sub-exposures to reject cosmic-ray events. The data were then 
wavelength and flux calibrated with conversion to heliocentric 
wavelengths and absolute flux units and rectified for the geo- 
metric distortions. The 2-D spectral image obtained for each 
slit position has a spatial scale of 0'.'0507 pixeL' along the slit, 
a dispersion of A/l = 1.108 A pixeL', and a spectral resolu- 
tion of H ^ 3000, covering the rest frame wavelength range 
6480-7050 A. 

For each spectrum we selected the regions containing 
the lines of interest. The lines were fitted, row by row, 
along the dispersion direction, together with a linear con- 
tinuum, with Gaussian functions using the task SPECFIT in 
STSDAS/IRAF. All emission lines present in the spectra (Ha, 
[N II]/1/16548,6583 and [S II]/l/t6716,6731) were fitted simulta- 
neously with the same velocity and width and with the relative 
flux of the [N II] lines kept fixed to 0.334. The assumption rela- 
tive to the sharing of the same kinematics for the emission lines 
will result in average kinematical quantities, weighted with the 
line surface brightnes|3- 

In the regions where the signal-to-noise ratio (SNR) was 
insufficient, the fitting was improved by co-adding two or more 
pixels along the slit direction. As an example of the quality of 
data, Fig.[T]shows the spectrum of NGC 5077, the only galaxy 
of the group with a marginal detection of a broad Ha compo- 
nent, at two different positions along NUC slit. 

3. Results 

The results obtained from the fitting procedure for the three 
galaxies are shown in Figs.|2]through|4l where we plot the kine- 
matical quantities, central velocity, and velocity dispersion, to- 
gether with the surface brightness of the emission fine with the 
best SNR at each location along the slit. The position-velocity 
diagram for the off-nuclear with the best quality of data is also 
reproduced in the upper panels of Figs. |2]and|4]for two sources 
of the program. In the following we discuss the results obtained 
for each galaxy of the sample. 

3.1. IC989 

The velocity curve on the nuclear slit shows a large-scale rota- 
tion with redshift at the negative space values (NE region) and 
blueshift at the positive (SW region) (Fig.|2]i. The same trend is 
observable in the off-nuclear slits, with almost the same values 
of velocity gradients. However, the velocity gradient is reversed 
at the very center, suggesting the presence of a counter-rotating 
nuclear gas disk. 

3.2. NGC 5077 

For the elliptical galaxy NGC 5077, the observed velocity 
curves are apparently consistent with gas in regular rotation 
around the galaxy's center. The velocity curve in the central 
slit, NUC (Fig.|3] central panel), has a full amplitude of ~ 400 
km s 'and shows a general reflection symmetry: starting from 



' However, we investigated the possibility that lines did not share 
common kinematics, but we did not find evidence of this. 
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Galaxy 


Type 


D 




N„i, 


PA 


Exp. T. 


Date 






(Mpc) 


(km s-') 






(s) 




IC 989 


E 


102 


185 


3 


-157 


494 


1998-03-13 


NGC 5077 


E3-4 


38 


256 


3 


-81 


418 


1998-03-12 


NGC 6500 


Sab 


42 


200 


5 


15 


236 


1998-11-03 



IC 989 



Row 601 





6700 

m 



Fig. 1. Nuclear spectra of NGC 5077: at the slit center position 
{upper panels), showing the possible presence of a weak broad 
Ha component, and at C'l from it {lower panels). The best 
model spectrum (dashed line) is superposed on the data; dotted 
lines are for the Gaussians and linear continuum components 
of the fit. Flux is in units of 10"'^ erg cm"^ s"' A"' arcsec"^. 
The residuals from the fits are shown in the smaller panels. 



the center, the velocity rises rapidly on both sides by ~ 200 
km s"' reaching a peak at r ~ (X'l from the center. This trend 
is followed by a small decrease and a substantially constant 
value at larger radii. Line emission and velocity dispersion are 
strongly peaked and smoothly decrease from the nucleus out- 
wards. 

The behavior seen in the off-nuclear slits is qualitatively 
similar to what is seen at the NUC location (Fig. |3] left and 
right panels), but with smaller velocity amplitude and, more 
important, a less extreme velocity gradient. Both the amplitude 
and gradient decrease at increasing distance of the slit center 
from the nucleus, with a behavior characteristic of gas rotating 
in a circumnuclear disk. 

Ground-based photomet r ic and k inematical studies 



dPemoulin-Ulrich et all 119841: iBertola et all Il991h show that 



NGC 5077 exhibits a gaseous disk with the major axis roughly 
orthogonal to the galaxy photometric major axis (PA ~ 10°). 
The gas isophotes show twisting and a marked warp on the 



> 




-0.4 -0.2 0.2 0.4 

Position along the slit [arcsecj 

Fig. 2. Upper panel: Velocity curve from the off-nuclear slit 
with the best quality of data. Lower panels: Velocity, surface 
brightness, and velocity dispersion along NUC slit. Surface 
brightness is in units of 10"'^ erg s"' 
along the slit are relative to the continuum peak. 



W side at r ~ 20"( ICaon et alj l2000l) . The gas has a fairly 
symmetric and smooth rotation curve at PA = 108° (~ major 
axis of the ionized gas distribution) with a half amplitude of 
~ 270 km s"' at r ~ 13". At PA = 10°, the stellar rotation curve 
exhibits a counter-rotating core {r < 5"). Along this axis, the 
gas rotates in the same direction as the stellar nucleus and 
shows a small-scale central velocity plateau. 



3.3. NGC 6500 

The velocity curves for this galaxy are complex (see Fig. |4|i. 
On the nuclear slit, the velocity apparently shows a general ro- 
tational trend, with redshift on the left and blueshift on the right 
side of the diagram. The same large-scale rotation is observable 
on the ofF-nuclear slit positions. Smaller scale high-amplitude 
(~ 200 kms"') velocity oscillations are also clearly visible. 

In the SW region at r < 0'.'4 the position-velocity diagrams 
of NUC and the adjacent SI slit show a trend that is inconsis- 
tent with gas in regular rotation around the galaxy's center (see 
Fig- in second upper panel). In this region the velocity curves 
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L.T 150 
^ 100 



Position along the slit [arcsec] 




Position along the slit [arcsoc 



Position along the slit [aresec] 



Fig. 3. Velocity, surface brightness, and velocity dispersion along the ofF-nuclear slit position Nl (left panel), NUC {central 
panel), and SI {right). Positions along the slits are relative to the continuum peak, positive values are NW. 



NGC 6500 



300 
200 



SI 



NUC 



Position along Lite sliL [arcsec] 

Fig. 4. Same as Fig|2] Second upper panel: Over-plot of NUC 
{empty squares) and SI {filled circles) velocity curves show- 
ing the inconsistency of the gas position-velocity diagram with 
regular rotation of a nuclear gas disk. 



least) a warped geometry for the disk, which cannot be con- 
strained with the present data. By these considerations, only 
NGC 5077 is a suitable candidate for providing a successful 
SMBH mass measurement. 



4. Modeling the rotation curves of NGC 5077 

described in detail in 



Our m odeling code, described in detail in iMarconi et al 
(l2003h . was used to fit the observed rotation curves of NGC 
5077. The code computes the rotation curves of the gas assum- 
ing that the gas is rotating in circular orbits within a thin disk in 
the galaxy potential. The gravitational potential has two com- 
ponents: the stellar potential (whose mass distribution will be 
determined in Sec. 14. Il l, characterized by its mass-to-light ratio 
and a dark mass concentration (the black hole), spatially unre- 
solved at HST+STIS resolution and characterized by its total 
mass Mbh- In computing the rotation curves, we take into ac- 
count the finite spatial resolution of the observations, and the 
intrinsic surface brightness distribution (ISBD) of the emission 
lines and we integrate over the slit and pixel area. The adopted 
HSTh-STIS point spread function (PSF) is obtained by fitting 
with three Gaussians the Tiny Tim (Krist & Hook 199 9) model 
PSF calculated at 6700 Awi thout the Lyot- stop, according to 
the procedure suggested by ' Pressell (|2006'). The is mini- 
mized to determin e the free paramet ers using the downhill sim- 
plex algorithm bv lPress et al.l (Il992h . 



diverge, reaching a separation of ~ 300 kms at r ~ 0'.'2. The 
over-plot of the position velocity diagrams suggests an expand- 
ing bubble of gas . 

Summarizing, we find that NGC 5077 shows a line emis- 
sion consistent with a circumnuclear gas disk in Keplerian ro- 
tation around the galaxy's center For the remaining two ob- 
jects, the trend of their velocity curves cannot be ascribed to 
the regular rotation of a nuclear gas disk. The position-velocity 
diagrams of IC 989 indicate the presence of a counter-rotating 
gas component. Modeling this configuration would require (at 



4.1. The Stellar mass distribution 

To assess the contribution of stars to the gravitational potential 
in the nuclear region, we derived the stellar luminosity den- 
sity from the observed surface brightness distribution. We re- 
constructed the galaxy light profile using a WFPC2 F547M {V 
band) image retrieved from the publ ic archive (Fig. [Sj. Based 
on WFPC2 F702W {R band) images. lTran et alj (2001) found a 
filamentary dust structure in NGC 5077, with position angle of 
the major axis of the main dust feature ~ 102° (at r - 10"). Our 
inspection of HST archive images, however, shows that the in- 
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Fig. 5. Left: Vband (F547M filter) image of NGC 5077. Right: Ha WFPC2/HST image. Field of view 4'.'5 x 4'.'5. North is up. 
East is left. 




Fig. 6. Results of the isophotes analysis of the V band image 
of NGC 5077. Surface brightness is shown in the top panel (in 
units of erg s"' cm"^ A"' arcsec"^), and the galaxy's ellipticity 
and position angle are shown in the middle and bottom panels, 
respectively. 



Fig. 7. Fit to the surface brightness profile obtained from an 
oblate spheroid stellar density distribution. 



Following Ivan der Marel & van den BoschI (Il998h . the stellar 
density distribution was parameterized as 



fluence of dust is negligible. In fact. iTran et al.l (120011) derived 
a mean visual extinction of only < Ay > - 0.048. We used the 
IRAF/STSDAS program ELLIPSE to fit elliptical isophotes to 
the galaxy (see Fig. |6]). Nuclear regions (r < 0'.'3) show the 
largest ellipticity variations. At larger radii, ellipticity and po- 
sition angle show small variations around the values 0.25 and 
10°, respectively. 

We derived the stars' density profile from the galaxy 
surface brightness following the same method described in 
Paper I when assuming an oblate spheroid density distribution. 



p(m) = po 



where m is given by m 



rb 



1 + 



rh 



■ + z?' Iq^, xyz is a reference 



system with the xy plane corresponding to the principal plane 
of the potential and q is the intrinsic axial ratio. The above 
stellar density is integrated along the line of sight, to provide 
the surface brightness distribution on the plane of the sky, and 
is then convolved with the telescope+instrument PSF (derived 
from Tiny Tim) and averaged over the detector pixel size to 
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obt ain the observed surf ace brightness distribution (for details 
see iMarconi et al.ll2003h . The de-projection was performed by 



adopting the reference value ofTy-l for the mass-to-light ra- 
tio in the V band (its real value will be determined in Sec. 14.2b . 
As the intrinsic potential axial ratio q is related to the inclina- 
tion of the principal plane with respect to the line of sight by 
the observed isophote ellipticity, we are left with the freedom 
of assuming different galaxy inclinations to the line of sight. 
We performed the de-projection adopting the value / - 44° for 
the galaxy inclination, as given from the HyperLeda database. 
Due to the small ellipticity of this galaxy, the precise value of 
its inclination only has a marginal effect on the resulting mass 
distribution. The best fit obtained is shown in Fig.|2]with a = 
0.47,yS = 0.79, and rb =0.59". 

4.2. Fitting the gas kinematics 

The observed kinematical quantities are averages over aper- 
tures defined by the slit width and the detector pixel size along 
the slit. The model fitting to the observed kinematical quantities 
thus depends on the intrinsic emission-line surface-brightness 
distribution, which is the weight for the averaging process, and 
on the following parameters: 

- the black hole mass Mbh; 

- the mass-to-light ratio of the stellar component, Tv; 

- the position of the kinematical center in the plane of the sky, 
xo, yo, with the reference frame centered on the nuclear slit 
and X axis aligned with the slit; 

- the inclination of the rotating disk, /; 

- the angle between the slits and the line of nodes, 0; 

- the systemic velocity of the disk, Vjys- 

With respect to Paper I, the modeling code has been slightly 
improved, so it now allows us to 1) fit the surface brightness 
and kinematical data simultaneously and 2) set limits to the 
parameters range. We used the second feature in particular to 
constrain t he ma ss-to-light ratio to Yy < 7.87, as derived from 
MarastonI (Il998h for a Salpeter initial mass function (IMF). 

A crucial issue for modeling the gas kinematics is the incli- 
nation of the nuclear gas disk, as this is strongly coupled with 
the black hole mass and Yy- In an oblate spheroid, the stable 
orbits of the gas are coplanar with the principal plane of the 
potential, and it is possible to directly associate the galaxy in- 
clination and line of nodes with those of the circumnuclear gas. 
However, the potential shape is not determined well enough 
by the isophotal fitting down to the innermost regions of the 
galaxy, and it is possible that a change of principal plane might 
occur at the smallest radii, in particular within the sphere of in- 
fluence of a supermassive black hole. We then performed a 
minimization for different values of /, namely / = 20°, ...,70°, 
allowing all other parameters to vary freely, including those de- 
scribing the intrinsic line surface brightness distribution. Due 
to the sensitivity of the velocity dispersions to the ISBD mod- 
eling and to other c omputational p roblems (i.e. a coarse sam- 
pling), as shown by Marconi et alJ 12006), we decided to ini- 
tially restrict the kinematical fitting only to the velocity curves. 
We show in Sec. l4.2.l1 that including the velocity dispersions in 
the fitting procedure only has a marginal effect on our results. 



Inspection of Fig. [3] suggests that the ISBD should be de- 
scribed at least by 2 components: one at the NUC slit (the 
nuclear component), plus a more extended one (base compo- 
nent). The intrinsic emission line surface brightness was mod- 
eled with a composition of two circularly symmetric Gaussians 
and a constant function (having as free parameters the ampli- 
tudes, scale radii, and the positions on the sky plane of the 
symmetry centers). iMarconi et al. ( 20061) have recently shown 
that the choice of a particular model for the line surface bright- 
ness distribution affects the final BH mass estimate very little, 
provided that the model reproduces the observed line emission 
within the errors. They showed that systematic errors on Mbh 
due to the adopted intrinsic emission line surface brightness 
were on the order of 0.05 in log Mbh for their STIS data. We 
also tested the effect, for our data, of using exponential func- 
tions in the ISBD modeling. The BH masses derived were $ 
0.03 in log Mbh lower than obtained using a Gaussian ISBD. 

From the best fits obtained at varying gas disk inclina- 
tion, we quoted the goodness of the fits only with respect 
to the velocity curves, since we are interested in the kine- 
matical model that reproduces the position-velocity data best. 
Excluding the contribution to the of the ISBD fit is a con- 
servative choice for our data, as we verified that quoting the 
global fits would narrow the ranges of acceptable model pa- 
rameters. The overall best-fitting model, presented in Fig. [8] is 
obtained for / = 40°, Mbh = 6.8 x IO^Mq, Yy = 7.7, and an 
angle from the slits to the line of nodes of 34° (i.e. offset by 25° 
from the ground-based measurements). Table|2]shows the best 
kinematical model parameters, together with the correspond- 
ing reduced value (x^ - ;^^/d.o.f. with d.o.f. - 16 degrees 
of freedom). Apparently, the best fit slightly underpredicts the 
large-scale velocity points. We then tested the effect of allow- 
ing a higher maximum value of the mass-to-light ratio, related 
to a different choice of the IMF. A fit was performed by adopt- 
ing the value 13.78 as limit to Yy, corresponding to the oldest 
stellar population and a giant-dominated IMF. Indeed, the best 
fit obtained (with Yy = 13.2) reproduces the observed veloci- 
ties better at large radii. However, the BH mass does not change 
significantly, resulting in log Mbh only 0.02 lower than our pre- 
vious result. We therefore preferred to maintain the constraint 
on Yy derived by adopting the standard Salpeter IMF. 

The value of minimum (see Table |2| is far higher than 
the value indicative of a good fit. To establish the statistical sig- 
nificance of v ariations among different models, we followed 
Barth et alJ (120011) and rescaled the error bars in our velocity 
measurements by adding a constant error in quadrature such 
that the overall best-fitting model provides ~ 1 . This is quite 
a conservative approach as it has the effect of increasing the fi- 
nal uncertainty on Mbh- The additional velocity error is 45 km 
s"'. The rescaledy^ values against inclination are reported in 
the bottom paneO of Fig.|9] The acceptable models (within the 
Icr confidence level, i.e. Ax^ < 0.06 for 16 d.o.f.), correspond 
to the range of inclination of 35° < / < 50°. 



- The different trend (high increase) of x^[scaled] at small angles 
with respect to NGC 3998 (Paper I, Fig. 11) is determined by the a 
priori fixing of an upper limit for Yy values. 
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Fig. 8. Overall best fit to the rotation curves (left) with the case of null BH mass (dashed line) superposed on the central slit. 
Corresponding fit to the line surface brightness distribution (right). From upper to bottom panel: Nl, NUC, SI. 



Table 2. Overall best-fit parameters. 
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The case of no black hole mass is also shown in Fig. 
[8] for the central slit. This model (obtained by constraining 
Ty < 7.87 and Vsys to the same value of the best fit) predicts 
a velocity gradient that is too shallow with respect to the data, 
particularly in the central region r < U.'l, with a correspond- 
ingly unacceptable value of - 7.9 (properly rescaled to the 
overall best fit). 

To evaluate the uncertainty associated to the black hole 
mass estimate, we explored its variation with respect to the pa- 
rameter that is more strongly coupled to it, i.e. the mass-to-light 
ratio Ty (having already considered the dependence on orienta- 
tion). The uncertainty on Mbh associated to changes in Ty has 
been estimated by building grid in the Mbh vs. Ty param- 
eter space. At each point of the grid, described by a fixed pair 
of Mbh and Tv values, we obtained the best fit model allowing 
all other parameters to vary freely and derived the correspond- 
ing Xr value (properly rescaled). The result of this analysis at 
i = 40° is presented in Fig. [10] The Icr range of the black hole 
mass at this disk inclination is Mbh = 6.^^^'^ x 10^ Mq, reported 
in the error bar in Fig. |9] 

If we conservatively adopt the same fractional uncertainties 
on Mbh, associated to variations in Ty, for the allowed disk 
inclinations, we obtain a global range of acceptable black hole 
mass of Mbh = (4.0 - 11.1) x lO^Mg at a Icr level. 



"X 1-2 



J . . I . . , tlL , . I . . . I J 

20 40 60 80 

i [degree] 

Fig. 9. Upper panel: Best fit values of Mbh obtained by vary- 
ing the disk inclination. The continuous line reproduces the oc 
1/sin^/ dependence of Mbh- Lower panel: Corresponding 
values scaled to the overall best fit. Inclinations smaller than 
35° or larger than 50° are excluded at a confidence level of Icr 
(shaded regions in the diagrams), evaluated for 16 degrees of 
freedom. 
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Fig. 10. contours at varying Ty and Mbh at an inclination 
of / = 40°. Contours are plotted for confidence levels of 1, 2, 
and 3cr. The plus sign marks the overall best fit. 



T 




-0.4 -0.2 0.2 0.4 

r [arcsec] 

Fig. 11. Velocity dispersion distribution expected from the 
overall best fit model (solid line) at ; - 40°, and derived 
(dashed line) by including the observed velocity dispersions in 
the fitting procedure. From upper to bottom panel: Nl, NUC, 
SI. 



4.2.1. Velocity dispersion distribution 

At this point in our analysis we tested the influence on the 
above results of including the velocity dispersions in the fit- 
ting procedure. The observed ionised-gas velocity dispersions 
are fairly large in NGC 5077, reaching almost 400 kms 'in 
the nuclear slit. Velocity dispersions larger than expected from 
unresolved rotation could be an indication of non circular mo- 



tions tha^£OuldJnvaM 



200ll:IVerdoes Kleiin et alj2000ll2002t[Cappellari et alj2002l) . 

However, in the case of NGC 5077, rotational and instrumental 
broadening is sufficient to reproduce the behaviour of velocity 
dispersions. 

First of all, by adopting the best fit parameter set (Table 
IJI, we obtained the velocity dispersion distribution shown in 
FigHH The observed velocity dispersions are reproduced ac- 
ceptably well by the model, with only a small mismatch of ob- 
served and model peak positions. Thus the nuclear rise of the 
velocity dispersion is accounted for as unresolved rotation by 
the fitting model. 

Furthermore, we repeated the minimization at / - 40°, 
this time by including the observed velocity dispersions in the 
fitting procedure (with no free parameters added). The best fit 
obtained is shown in Fig. [TT] An improvement of the match 
with observed values is obtained, while there is a slight wors- 
ening of the fit to the velocity curves, with Xi (scaled to the 
overall best fit) - 1.08. The effect of including the line widths 
in the analysis results in increasing the BH mass by only ~ 24% 
to Mbh = 8.4 x IO'^Mq. 



5. Discussion 

The model fitting performed for NGC 5077 indicates that 
the kinematics of gas in its innermost regions is described 
correctly as circular motions in a thin disk when a point- 
like dark mass of Mbh = 6.8+^-^ x IO^Mq is added to 
the galaxy potential. We can now explore how BH mass re- 
lates to the properties of the host ga laxy, including spheroid 
(bulg e) massdM arconi & junti 2003)). stellar veloc ity disper- 
sion (iTremaine et al.ll2002t iFerrarese & Fordll2005h . and light 
concentration (Graham & Driver 2007 ). 

Following lMarconi & Hunt (l2003h. we used the virial mass 
Mvir = kReO-ljG, with k ^ 5 ("Capp ellari et al.ll200^ to de- 
termine the bulge mass of NGC 5077. We found discrepant 
values for the effective radius of the bulge. Re, in the litera- 
ture. The values were derived through different methods (i.e. 
isophotal fitting, photometric aperture growth cu rves) and ar e 



15'.'4 (Bender etal. 19891), 15'.'9 JSanchez-Portal et al 



17.'3 (Poulain&Nieto 1994), 19"(Bertola et al. 1991), 22'.'1 
(iTruiillo et al...2004). 22V8 (de Vaucouleurs et al...l99ll RC3), 



i 



2004) 



and 25"( lBurstein et al.lll987l) . We decided to adopt the average 
of the above values, R^. - 19'.'6 + 3'.'7 as the effective bulge ra- 
dius of NGC 5077. At the adopted distance of the galaxy, this 
value corresponds to 3.6 + 0.7 kpc. 

The stellar velocity dispersion d eterminations we fo und in 
the literature are 285 



28 kms '([Davies et al. 
4"aperture), 252 + 12 km s ' jCarollo etal.lll993 



1987, 1'.'5 X 



1'.'3 X 5"and 



1'.'6 X 5"apertures), and 263 + 5 km s ' fCaon et allbOOOi 1'.'5 
X r.'8 aperture). All the above values refer to the very cen- 
tral regions of the galaxy, corresponding to equivalent circu- 
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Fig. 12. Mbh vs bulge mass from lMarconi & Hunt! (l2003h with 
the best fit obtained from a bisector linear regression analysis 
(soHd line). The positions of NGC 5077 and NGC 3998 are 
shown by empty circles. 



lar apertures (see ICappellari et al. I l2006h from ~ l"to I'.'S. We 
then normahzed the values of the central velocity dispersion to 
a circular aperture of radius equal to the adopted Rg, following 
the formula for aperture corrections by Cappellari et al.. (2006) . 
The average of these normalized values leads to ctr^ = 222 +15 
kms~'. 

Using this value of cr as an approximation for cTbui, we ob- 
tained the value Mbui - (2.1 +0.7)x 10"Mp for t he bulge mass. 
From the correlation of Marconi & Hunt (12003 ^. the expected 
Mbh for NGC 5077 is 3.0 xIO^Mq, in agreement within a fac- 
tor ~ 2.3, with our determination (see FigfTSli. 

For the correlations between BH mass and central stellar 
velocity dispersion, o ne obtains 2.0 xlO^M g for NGC 5077 
adopting the form by iTremaine et alJ (l2002h (see Fig[T3] up- 
per panel), a factor 3.4 lower than our estimate. We also 
compared our measurement with the expected value from the 
Mbh - cr relation derived by Ferrarese & Ford (2005 E Mbh = 
5.4 xIO^Mq, in full agreement with our estimate (see FigfT3l 
bottom panel). 

Conversely, the BH mass expected for NGC 5077 on the 
basis of its correlation with the Sersic concentration index 
(from .Trujillo et al.. 2004) is a factor 5~6 lower than our es- 
timate when using the l inear Mbh- Sersic index correlation 
(IGraham & Driver hoof ) and a factor 3 lower when using the 
qu adratic form they propos ed. 



Marconi & Hunt! (120031) noticed that Mbh is separately cor- 



related with both cr and Rg. This is shown by the weak corre- 
lation between the residuals of the Mbh - cr correlation with 



^ We normalized the central velocity dispersion to an aperture of 
radius equal to 1/8 of and derived cr^/g = 255 ± 17 km s"' . 



10 



9 - 



^8 - 



If 

G 



O 



7 - 



1 



1.8 2 2.2 2.4 2.6 
Log a [km s"'] 



10 



SB - 



o 




Fig. 13. Mbh vs central 
from Tremaine et al. 



i 



errarese 



&Fordl (l2005l) (/, 



stellar velocity dispersion 
(120021) (upper panel) and from 
'ower panel). Solid lines indi- 



cate the best fits for the two correlations. The positions of 
NGC 5077 and NGC 3998 are shown by empty circles. 



Re reproduced here in Fig. [14] for the Tremaine at al. fornfl 
Our determinations of the black hole mass in NGC 5077 and 
in NGC 3998 (Paper I) support this idea. Unlike NGC 3998, 
which has one of the lowest values of 7?e among galaxies with 
measured Mbh (0.85 kpc) and shows a negative residual from 
the Mbh - o" correlation, NGC 5077 has an intermediate R^ 
value (3.6 kpc) and show a small but positive residual (see Fig. 
fl4l i. In the same sense, but with a larger pos itive residual, there 
is the result found by ICapetti et alT(l2005l) for NGC 5252, a 
galaxy with quite a large effective radius (9.7 kpc). This indi- 



A similar trend is obtained using the correlation in the Ferrarese 
et al. form. 
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Fig. 14. Residuals from the Mbh vs cr correlation (in the 
Tremaine et al. form) reported against the bulge effective radius 
^e- The positions of NGC 5077 and NGC 3998 are indicated by 
empty circles. 

cates the presence of a black hole's "fundamental plane" in the 
sense that a combination of at least cr and drives the cor- 
relations between Mbh and the bulge properties. The physical 
implications of these results will be discussed in a forthcoming 
paper, when new direct BH mass measurements are presented, 
in order to base our discussion on a statistically more significa- 
tive sample of Mbh determinations. 

6. Summary and conclusions 

We have presented results from a gas kinematics study in the 
nucleus of three neai'by LINER galaxies: IC 989, NGC 5077, 
and NGC 6500 using archival HST/STIS spectra. Only in the 
case of NGC 5077 does the nuclear velocity curves appear to 
be associated with gas in regular rotation. For IC 989 our re- 
sults indicate an inner counter-rotating gas system, while for 
NGC 6500 the complex trend in the velocity curves suggests a 
nuclear expanding bubble. 

We used our modeling code to fit the observed [N II]/16583 
surface brightness distribution and velocity curve of NGC 
5077. The dynamics of the rotating gas can be accurately re- 
produced by motions in a thin disk when a compact dark mass 
of Mbh = 6.8^28 ^ 10** Mq is added to the stellar mass com- 
ponent. Furthermore, the black hole in NGC 5077 has a sphere 
of influence radius, /?sph - GM^iijcrl^.^^, of ~ 62 pc 0'.'34), 
well-resolved at the HST resolution (2/?sph/^res - 6.8). 

For what concerns the connections of this BH mass esti- 
mate with the properties of the host galaxy, the Mbh value for 
NGC 5077 is in good agreement (within a factor of 2.3) with 
the Mbh - A^bui correlation between BH and host bulge mass. 
The black hole mass predicted by the Mbh - cr correlation is 



a factor of 3. 4 lower than our measure, adopting the relation 
found by FTremaine et alJ (l2002h and in excellent agre ement us- 
ing the parameterization by 'F errarese & FordI ( 2005 ). 

This result, in conjunction with the previous results for 
NGC 3998 (Paper I), strengthens the possibility of a connection 
between the residuals from the Mbh - cr relation and the bulge 
effective radius. While NGC 3998, indeed, has one of the low- 
est values of among galaxies with measured Mbh and shows 
a negative residual, NGC 5077 has a larger effective radius and 
shows a small positive residual. We also recently showed that 
the same result was found for the Seyfert galaxy NGC 5252: a 
larger effective radius corresponds in this galaxy to a still larger 
positive residual. 

Apparently, only with a combination of at least cr and 
is it possible to account for the correlations between Mbh and 
other bulge properties, indicating the presence of a black hole's 
"fundamental plane". Clearly, the number of direct black-hole 
mass measurements must be further increased, together with 
precise determinations of cr and R^, to test these conclusions 
on a stronger statistical basis. In a forthcoming paper we wiU 
present new BH mass determinations and discuss the physical 
implications of our results. 
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